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We consider the impact of imposing generalized CP symmetries on the Higgs sector of the two- 
Higgs-doublet model, and identify three classes of symmetries. Two of these classes constrain the 
scalar potential parameters to an exceptional region of parameter space which respects either a Z2 
discrete flavor symmetry or a (7(1) symmetry. We exhibit a basis-invariant quantity that distin- 
guishes between these two possible symmetries. We also show that the consequences of imposing 
these two classes of CP symmetry can be achieved by combining Higgs Family symmetries, and 
' (-H ' , that this is not possible for the usual CP symmetry. We comment on the vacuum structure and on 

^ • renormalization in the presence of these symmetries. Finally, we demonstrate that the standard CP 

I ' symmetry can be used to build all the models we identify, including those based on Higgs Family 

symmetries. 

PACS numbers: ll.30.Er, 12.60.Fr, 14.80.Cp, ll.30.Ly 

> ■ I. INTRODUCTION 

l> 

' Despite the fantastic successes of the Standard Model (SM) of electroweak interactions, its scalar sector remains 
largely untested An alternative to the single Higgs doublet of the SM is provided by the two-Higgs-doublet 
• ' model (THDM), which can be supplemented by symmetry requirements on the Higgs fields $1 and $2- Symmetries 
, leaving the kinetic terms unchanged^ may be of two types. On the one hand, one may relate <I>a with some unitary 
transformation of <i>t,. These are known as Higgs Family symmetries, or HF symmetries. On the other hand, one 
\ may relate $a with some unitary transformation of <i>^. These are known as generalized CP symmetries, or GCP 
• • symmetries. In this article we consider all such symmetries that are possible in the THDM, according to their impact 
on the Higgs potential. We identify three classes of GCP symmetries. 

The study is complicated by the fact that one may perform a basis transformation on the Higgs fields, thus hiding 
\^ ■ what might otherwise be an easily identifiable symmetry. The need to seek basis invariant observables in models 
5^ , with many Higgs was pointed out by Lavoura and Silva Q, and by Botella and Silva @, stressing applications to 
CP violation. Refs. 0,0] indicate how to construct basis invariant quantities in a systematic fashion for any model, 
including multi-Higgs-doublet models. Work on basis invariance in the THDM was much expanded upon by Davidson 
and Haber [8|, by Gunion and Haber by Haber and O'Neil [iiIj and by other authors T he p revious 

approaches highlight the role played by the Higgs fields. An alternative approach, spearheaded by Nishi by 
Ivanov [1, 3 and by Maniatis et al ITsI . hig hlights the role played by field bilinears, which is very useful for studies of 
the vacuum structure of the model jl6l. [l7|. In this paper, we describe all classes of HF and GCP symmetries in both 
languages. One problem with two classes of GCP identified here is that they lead to an exceptional region ofparameter 
space (ERPS) previously identified as problematic by Gunion and Haber ^ and by Davidson and Haber [8|. Indeed, 



X 



^ It has been argued by Ginsburg Q and by Ivanov 0, U that one should also consider the effect of non-unitary global symmetry 
transformations of the two Higgs fields, as the most general renormalizable Higgs Lagrangian allows for kinetic mixing of the two Higgs 
fields. In this work, we study the possible global symmetries of the efl'ective low-energy Higgs theory that arise after diagonalization 
of the Higgs kinetic energy terms. The non-unitary transformations that diagonalize the Higgs kinetic mixing terms also transform the 
parameters of the Higgs potential, and thus can determine the structure of the remnant Higgs flavor symmetries of effective low-energy 
Higgs scalar potential. It is the latter that constitutes the main focus of this work. 
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no basis invariant quantity exists in the literature that distinguishes between the Zi and f/(l) HF symmetries in the 
ERPS. 

If evidence for THDM physics is revealed in future experiments, then it will be critical to employ analysis techniques 
that are free from model-dependent assumptions. It is for this reason that a basis-independent formalism for the 
THDM is so powerful. Nevertheless, current experimental data already impose significant constraints on the most 
general THDM. In particular, we know that custodial symmetry breaking effects, flavor changing neutral current 
(FCNC) constraints, and (to a lesser extent) CP-violating phenomena impose some significant restrictions on the 
structure of the THDM (including the Higgs-fermion interactions). For example, the observed suppression of FCNCs 
implies that either the two heaviest neutral Higgs bosons of the THDM have masses above 1 TeV, or certain Higgs- 
fermion Yukawa couplings must be absent [Tsl . The latter can be achieved by imposing certain discrete symmetries 
on the THDM. Likewise, in the most general THDM, mass splittings between charged and neutral Higgs bosons can 
yield custodial-symmetry breaking effects at one-loop that could be large enough to be in conflict with the precision 
electroweak data jl9j . Once again, symmetries can be imposed on the THDM to alleviate any potential disagreement 
with data. The implications of such symmetries for THDM phenomenology has recently been explored by Gerard and 
collaborators [20^ and by Haber and O'Neil [2l'|. 

Thus, if THDM physics is discovered, it will be important to develop experimental methods that can reveal the 
presence or absence of underlying symmetries of the most general THDM. This requires two essential pieces of input. 
First, one must identify all possible Higgs symmetries of interest. Second, one must relate these symmetries to basis- 
independent observables that can be probed by experiment. In this paper, we primarily address the first step, although 
we also provide basis-independent characterizations of these symmetries. Our analysis focuses the symmetries of the 
THDM scalar potential. In principle, one can extend our study of these symmetries to the Higgs-fermion Yukawa 
interactions, although this lies beyond the scope of the present work. 

This paper is organized as follows. In section |TT] we introduce our notation and define an invariant that does 
distinguish the Z2 and C/(l) HF symmetries in the ERPS. In section UIIl we explain the role played by the vacuum 
expectation values in preserving or breaking the C/(l) symmetry, and we comment briefly on renormalization. In 
section lTVl we introduce the GCP transformations and explain why they are organized into three classes. We summarize 
our results and set them in the context of the existing literature in section |Vl and in section IVTl we prove a surprising 
result: multiple applications of the standard CP symmetry can be used to build all the models we identify, including 
those based on HF symmetries. We draw our conclusions in section IVlIl 



II. THE SCALAR SECTOR OF THE THDM 



A. Three common notations for the scalar potential 



Let us consider a S'C/(2) ® C/(l) gauge theory with two Higgs-doublets <f>a, with the same hypercharge 1/2, and 
with vacuum expectation values (vevs) 





fa/%/2 



(1) 



The index a runs from 1 to 2, and we use the standard definition for the electric charge, whereby the upper components 
of the SU{2) doublets are charged and the lower components neutral. 
The scalar potential may be written as 

Vh = mii^l^i + m22^'2$2 - mi2'I'I$2 + H.c. 

+ lAi($t$l)2 + iA2($t$2)' + A3($l$i)($5$2) + A4($I$2)($|$l) 

+ [5A5($1$2)' + A6($I$l)($I$2) + A7($J$2)($1*2) + H.C.] , (2) 

where m^ii "^22 j ^^'^ Ai, • • • , A4 are real parameters. In general, m^2i -^Si '^6 and A7 are complex. "H.c." stands for 
Hermitian conjugation. 

An alternative notation, useful for the construction of invariants and championed by Botella and Silva Q is 



Vh = rafc($l$fc) + iZ,b,c<i($I*6)($I*d), 



(3) 
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where Hermiticity implies 



b.cd 



■^cd.ab 



V 



ha 1 



= z, 



ha Ac' 



The extremum conditions are 



[Yab + Zab,cd VdVc] Wfc = 



(for a = 1,2) 



Multiplying by v* leads to 



Yab{vlVb) = -Zab.cd (vlvb) (w>c)- 



(4) 
(5) 
(6) 



One should be very careful when comparing Eqs. ^ and ^ among different authors, since the same symbol may 
be used for quantities which differ by signs, factors of two, or complex conjugation. Here we follow the definitions of 
Davidson and Haber M. With these definitions: 



and 





Yn 


= mil, 


Y12 = 








Y21=- 




Y22 = 


- '^22 > 






Z\\.\\ 


= Ai, 


^22,22 


= A2, 




■^11,22 


= -^22,11 


= A3, 


^12,21 


= ^21,12 


— A4 




^12,12 


= As, 


^21,21 


- ■^5- 




■^11,12 


= ^12,11 


— Ae, 


^11,21 


= ^21,11 


-a; 


■^22,12 


= Z\2.22 


= A7, 


^22,21 


= ^21,22 


= a; 



(7) 



(8) 

The previous two notations look at the Higgs fields $a individually. A third notation is used by Nishi [l^, [ij] and 
Ivanov 01 , who emphasize the presence of field bilinears {^l^b) [17]. Following Nishi 1^ we write: 



(9) 



where /i = 0, 1, 2, 3 and 





1 
2 




($t$2) 
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1 
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($I$2) + 


($t$i) 




^2 = 


I 
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($1$2) 






1 
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"(<i>I<l'i) - 


($t$2) 





(10) 

In Eq. ([5]), summation of repeated indices is adopted with Euclidean metric. This differs from Ivanov's notation [1,0], 
who pointed out that parametrizes the gauge orbits of the Higgs fields, in a space equipped with a Minkowski 
metric. 

In terms of the parameters of Eq. ([2]), the 4- vector and 4x4 matrix A^t^ are written respectively as: 



'-22' 



'2Rem^2: 2Im 



'12' 



^22 ) ; 



and 



Ku = 



(Ai+A2)/2 + A3 Re(A6 + A7) -Im(A6 + A7) 
Rc(A6 + A7) A4 + ReAs — ImAs 

— Im(A6 + A7) —ImAs A4 — ReAs 

(Ai-A2)/2 Re(A6-A7) -Im(A6-A7) 



(Ai - A2)/2 
Rc (Ag - A7) 
-Im(A6 - A7) 
(Ai+A2)/2-A3 



Eq. (O is related to Eq. ^ through 



A^" = ^Zab.cdCTba^dc^ 

where the matrices cr* are the three Pauli matrices, and cr*^ is the 2x2 identity matrix. 



(11) 



(12) 



(13) 
(14) 
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B. Basis transformations 

We may rewrite the potential in terms of new fields obtained from the original ones by a simple (global) basis 
transformation 

^a^K = Uab'i'b, (15) 

where U £ U{2) is a 2 x 2 unitary matrix. Under this unitary basis transformation, the gauge-kinetic terms are 
unchanged, but the coefficients Yat and Zat.cd are transformed as 

Yab -> Y^,^UaaY^0U:i,, (16) 
Zab.cd ^'ab.cd = UaaUcj ZafS.jsU^i^U^g, (17) 



and the vevs are transformed as 



Va ^ v'^ = UabVb- (18) 



Thus, the basis transformations U may be utilized in order to absorb some of the degrees of freedom of Y and/or Z, 
which implies that not all parameters of Eq. ([Sj have physical significance. 

C. Higgs Family symmetries 

Let us assume that the scalar potential in Eq. ([3]) has some explicit internal symmetry. That is, we assume that 
the coefficients of Vh stay exactly the same under a transformation 

$a ^ $f = Sab^b- (19) 

5 is a unitary matrix, so that the gauge-kinetic couplings are also left invariant by this Higgs Family symmetry (HP 
symmetry). As a result of this symmetry, 

Yab — Yi^Ij — Saa Y^p S/^p, (20) 
Zab,cd = Z^j^^^j^ = Saa Sc-f Zaf}^^S Sip S*i^. (21) 

Notice that this is not the situation considered in Eqs. p ^ -(fT7 |) . There, the coefficients of the Lagrangian do 
change (although the quantities that are physically measurable are invariant with respect to any change of basis). In 
contrast, Eqs. p9|) - ((2T|) imply the existence of a HF symmetry S of the scalar potential that leaves the coefficients of 
Vh unchanged. 

The Higgs Family symmetry group must be a subgroup of full U (2) transformation group of 2 x 2 unitary matrices 
employed in Eq. p^ . Given the most general THDM scalar potential, there is always a U(l) subgroup of U{2) under 
which the scalar potential is invariant. This is the global hypercharge U{1)y symmetry group: 

U{1)y: $i^e*''$i, $2^e^^$2, (22) 

where 9 is an arbitrary angle (mod 27r). The invariancc under the global U (l)y is trivially guaranteed by the invariance 
under the SU{2) (g) C/(l) clcctrowcak gauge symmetry. Since the global hypercharge U{1)y is always present, we shall 
henceforth define the HF symmetries as those Higgs Family symmetries that are orthogonal to U{1)y. 

We now turn to the interplay between HF symmetries and basis transformations. Let us imagine that, when written 
in the basis of fields $a, Vh has a symmetry S. We then perform a basis transformation from the basis ^a to the basis 
as given by Eq. (|15p . Clearly, when written in the new basis, Vh does not remain invariant under S. Rather, it 
will be invariant under 

S' = USUI (23) 

As we change basis, the form of the potential changes in a way that may obscure the presence of a HF symmetry. In 
particular, two HF symmetries that naively look distinct will actually yield precisely the same physical predictions if 
a unitary matrix U exists such that Eq. (I23|) is satisfied. 

HF symmetries in the two-Higgs-doublet model (THDM) have a long history. In papers by Glashow and Weinberg 
and by Paschos [TSj], the discrete Z2 symmetry was introduced, 

Z2: $1 ^ $1, $2 ^ -$2, (24) 
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in order to preclude flavour-changing neutral currents This is just the interchange 

Ha : $1 ^ $2, (25) 
seen in a different basis, as shown by applying Eq. ()23|) in the form 

Peccei and Quinn [22( | introduced the continuous C/(l) symmetry 

U{1): $i^e-*«$i, $2^e'^$2, (27) 

true for any value of 9, in connection with the strong CP problem. Of course, a potential invariant under U{1) is also 
invariant under Z2. 

Finally, we examine the largest possible Higgs Family symmetry group of the THDM, namely U{2). In this case, 
a basis transformation would have no effect on the Higgs potential parameters. Since Sab is the only t/(2)-invariant 
tensor, it follows that 

Yab - CiSab, (28) 
Zab,cd = C26abScd + CzSad^bc , (29) 



where ci, C2 and C3 are arbitrary real numbers. ^ One can easily check from Eqs. (jl6p and (|17p that the unitarity of 
U implies that Y' = Y and Z' ^ Z for any choice of basis, as required by the ?7(2)-invariance of the scalar potential. 
Eqs. (pS)) and (P^jl impose the following constraints on the parameters of the THDM scalar potential (independently 
of the choice of basis): 

22 2 n 

77122 = "111, "^12 = 0: 

Ai = A2 = A3 + A4 , As = Ae = Ay = . (30) 

As there are no non-zero potentially complex scalar potential parameters, the [/(2)-invariant THDM is clearly CP- 
invariant. 

As previously noted, the U{2) symmetry contains the global hypercharge U{1)y as a subgroup. Thus, in order to 
identify the corresponding HE symmetry that is orthogonal to C/(l)y, we first observe that 

U{2) = SU{2) (g> U{1)y/Z2 = SO{S) ® U{1)y ■ (31) 

To prove the above isomorphism, simply note that any U{2) matrix can be written as C/ = e'^C/, where U S SU{2). 
To cover the full U{1)y group, we must take < 9 < 2tt. But since both U and —U are elements of SU{2) whereas 
-1-1 and — 1 = e^^ are elements of U{\)y, we must identify U and —U as the same group element in order not to 
double cover the full U{2) group. The identification of U with —U in SU{2) is isomorphic to 5*0(3), using the well 
known isomorphism 5*0(3) = SU{2)/Z2. Consequently, we have identified S0(3) as the HE symmetry that constrains 
the scalar potential parameters as indicated in Eg. ([50)1 . 

The impact of these symmetries on the potential parameters in Eq. ^ is shown in section|Vl As mentioned above, 
if one makes a basis change, the potential parameters change and so does the explicit form of the symmetry and of 
its implications. For example, Eq. shows that the symmetries Z2 and H2 are related by a basis change. However, 
they have a different impact on the parameters in their respective basis. This can be seen explicitly in Table U of 
section fVl One can also easily prove that the existence of either the Z2, H2 or Peccei-Quinn U{1) symmetry is 
sufficient to guarantee the existence of a basis choice in which all scalar potential parameters are real. That is, the 
corresponding scalar Higgs sectors are explicitly CP-conserving. 

Basis invariant signs of HE symmetries were discussed extensively in Ref. Q. Recently, Ferreira and Silva [2^ 
extended these methods to include Higgs models with more than two Higgs doublets. 

Consider first the THDM scalar potentials that are invariant under the so-called simple HE symmetries of Ref. [1^ . 
We define a simple HE symmetry to be a symmetry group G with the following property: the requirement that the 
THDM scalar potential is invariant under a particular element g (z G (where g ^ e and e is the identity element) is 



^ Note that there is no Sac^bd term contributing to Z^h.cd^ such a term is not invariant under the transformation of Eq. II17I I. 
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sufRcient to guarantee invariance under the entire group G. The discrete cycUc group Zn = {e , g , , ■ ■ ■ , 
where = e, is an example of a possible simple HF symmetry group. If we restrict the TDHM scalar potential to 
include terms of dimension- four or less (e.g., the tree-level scalar potential of the THDM), then one can show that 
the Peccei-Quinn U{1) symmetry is also a simple HF symmetry. For example, consider the matrix 

^ = ( e^W3 ) . (32) 

Note that S is an element of the cyclic sub-group Z^^ — {S , S'^ , — 1} of the Peccei-Quinn U{1) group. As shown 
in Ref. [1^, the invariance of the tree-level THDM scalar potential under $a Sab^b automatically implies the 
invariance of the scalar potential under the full Peccei-Quinn U{1) group. In contrast, the maximal HF symmetry, 
50(3), introduced above is not a simple HF symmetry, as there is no single element of S* G 50(3) such that 
invariance under <l>a Sab^b guarantees invariance of the tree- level THDM scalar potential under the full SO (3) 
group of transformations. 

Typically, the simple HF symmetries take on a simple form for a particular choice of basis for the Higgs fields. We 
summarize here a few of the results of Ref. jl^l : 

1. In the THDM, there are only two independent classes of simple symmetries: a discrete Z2 flavor symmetry, and 
a continuous Peccei-Quinn U{1) flavor symmetry. 

2. Other discrete flavor symmetry groups G that are subgroups of U{1) are not considered independent. That is, 
if 5 € G (where 5* 7^ e), then invariance under the the discrete symmetry $ — s- 5$ makes the scalar potential 
automatically invariant under the full Peccei-Quinn C/(l) group; 

3. In most regions of parameter space, one can build quantities invariant under basis transformations that detect 
these symmetries; 

4. There exists a so-called exceptional region of parameter space (ERPS) characterized by 

22 2 A 

A2 = Ai, Ay = -Ag. (33) 

As shown by Davidson and Haber Q , a theory obeying these constraints does have a Z2 symmetry, but it may 
or not have a U{1) symmetry. Within the ERPS, the invariants in the literature cannot be used to distinguish 
the two cases. 

The last statement above is a result of the following considerations. In order to distinguish between Z2 and f/(l), 
Davidson and Haber ^ construct two invariant quantities given by Eqs. (46) and (50) of Ref. [8j]. Outside the ERPS, 
these quantities are zero if and only if U{1) holds. Unfortunately, in the ERPS these quantities vanish automatically 
independently of whether or not U{1) holds. Similarly, Ferreira and Silva psj have constructed invariants detecting 
HF symmetries. But their use requires the existence of a matrix, obtained by combining Yab and Zab,cd, that has two 
distinct eigenvalues. This does not occur when the ERPS is due to a symmetry. Finally, in the ERPS, Ivanov Q 
states that the symmetry might be "(^2)^ or 0(2)" [our Z2 or our C/(l)] and does not provide a way to distinguish 
the two possible flavor symmetries [2J|. 

Gunion and Haber 9] have shown that the ERPS conditions of Eq. p3|) are basis independent; if they hold in one 
basis, then they hold in any basis. Moreover, for a model in the ERPS, a basis may be chosen such that all parameters 
are real.'^ Having achieved such a basis, Davidson and Haber Q demonstrate that one may make one additional basis 
transformation such that 

22 2 n 

™22=™ii> m-12 = 0) 

A2 = Ai, A7 = A6=0, ImA5 = 0. (34) 

These conditions express the ERPS for a specific basis choice. 

One might think that this is such a special region of parameter space that it lacks any relevance. However, the fact 
that the conditions in Eq. p3p hold in any basis is a good indication that a symmetry may lie behind this condition. 



^ Given a scalar potential whose parameters satisfy the ERPS conditions with Im(A5Ag) ^ 0, the unitary matrix required to transform 
into a basis in which all the scalar potential parameters are real can be determined only by numerical means. 
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Indeed, as pointed out by Davidson and Haber [8|, combining the two symmetries Z2 and 112 same basis one 

is lead immediately to the ERPS in the basis of Eq. Up to now, we considered the impact of imposing on the 

Higgs potential only one symmetry. This was dubbed a simple symmetry. Now we are considering the possibility 
that the potential must remain invariant under one symmetry and also under a second symmetry; this implies further 
constraints on the parameters of the Higgs potential. We refer to this possibility as a multiple symmetry. As seen 
from Table HI of section fVl imposing Z2 and 112 in the same basis leads to the conditions in Eq. (l34|) . Incidentally, this 
example shows that a model which lies in the ERPS, is automatically invariant under Z2. 

In section ITVl we will show that all classes of non-trivial CP transformations lead directly to the ERPS, reinforcing 
the importance of this particular region of parameter space. 



D. Requirements for U{1) invariance 

In the basis in which the U{1) symmetry takes the form of Eq. ([27ll . the coefficients of the potential must obey 

m[l = 0, = = A'7 = 0. (35) 

Imagine that we have a potential of Eq. ([2]) in the ERPS: mfj^ = '^22^ "^12 = 0; A2 = Ai, and A7 = We now 

wish to know whether a transformation U may be chosen such that the potential coefficients in the new basis satisfy 
the U{1) conditions in Eq. (|35p. Using the transformation rules in Eqs. (A13)-(A23) of Davidson and Haber Q, we 
find that such a choice of U is possible if and only if the coefficients in the original basis satisfy 

2A3-A5A6(Ai-A3-A4)-A2a;-0, (36) 

subject to the condition that AgAg is real. 



E. The D invariant 



Having established the importance of the ERPS (as it can arise from a symmetry), we will now build a basis 
invariant quantity that can be used to detect the presence of a U(l) symmetry in this special case. 

The quadratic terms of the Higgs potential are always insensitive to the difference between Z2 and U{1). Moreover, 
the matrix Y is proportional to the unit matrix in the ERPS. One must thus look at the quartic terms. We were 
inspired by the expression of A^i^ in Eq. (fT2|) . which appears in the works of Nishi [l^ [3| and Ivanov [1, In 
the ERPS of Eq. ([33]), A^^ breaks into a 1 x 1 block (Aqo), and a 3 x 3 block (A = {A^}; i,j = 1,2,3). A basis 
transformation U belonging to SU{2) on the fields corresponds to an orthogonal 5*0(3) transformation in the 
bilinears, given by 

Oij = ^TT[U^a,Uc7,]. (37) 

Any matrix O of S0{3) can be obtained by considering an appropriate matrix U of SU{2) (unfortunately this 
property does not generalize for models with more than two Higgs doublets). A suitable choice of O can be made 
that diagonalizes the 3x3 matrix A, thus explaining Eq. ([M]) . In this basis, the difference between the usual choices 
for C/(l) and Z2 corresponds to the possibility that ReAs might vanish or not, respectively. 

We will now show that, once in the ERPS, the condition for the existence of U{1) is that A has two eigenvalues 
which are equal. The eigenvalues of a 3 x 3 matrix are the solutions to the secular equation 

^3 



a2X 



aix + oo = 0, 



(38) 



where 



ao 



iTr (A^) - i(Tr A)3 + i(Tr A)Tr(A2) 

^Zab,cdZdc,ba'^1^ 



z 



dc,gh 
3 



z, 



hg.ba 



2 

_ 37(2)7(2) 
2^dc ^ba 



ai 



02 



det A = - 

_ 1 7 

3 ^ab.cd 

TrZ«V + ifTrZ(i 

i(TrA)2 - iTr(A2) 
\ 2 

TrZ^^M -TrZ(i'TrZ(2 



ZH) _ 1^(2) 



Tr 



b.cd^dc.ba 



-TrA 



(39) 



(40) 



(41) 
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and 



7(1) = 7 _ / Ai + A4 Ae + A7 \ , . 

^ab - ^aa.ah - ^ XI + X*j A2 + A4 / ' ^ ^ 

7(2) _ 7 _ Ai + A3 Ae + A7 A / .„x 

= ^aa.ab - I + A? A2 + A3 j ■ ^^"^^ 



The cubic equation, Eq. has at least two degenerate solutions if [1^ 

^ = [iai - lalf+ [i(aia2 - 3ao) - 2y«2]' (44) 

vanishes. 

The expression of D in terms of the parameters in Eq. ([2]) is rather complicated, even in the ERPS. But one can 
show by direct computation that if the C/(l)-symmetry condition of Eq. holds (subject to AgApbeing real), then 
D = 0. We can simplify the expression for D by changing to a basis where all parameters are real [9|, where we get 

D = -^ [A5(Ai - A3 - A4 + As) - 2X1] ' [(Ai - A3 - A4 - As)^ + IGA^] . (45) 
If Ag 7^ 0, then D = means 

2A2 = A5(Ai-A3-A4 + A5). (46) 

If Ae = 0, then £> = corresponds to one of three possible conditions: 

As^O, A5 = ±(Ai-A3-A4). (47) 



Notice that Eqs. (|46|) and ([47]) are equivalent to Eq. ((36)) in any basis where the coefficients are real. 

Although D can be defined outside the ERPS, the condition D = only guarantees that the model is invariant 
under U{1) inside the ERPS of Eq. I|33p. Outside this region one can detect the presence of a U{1) symmetry with the 
invariants proposed by Davidson and Haber Q . This closes the last breach in the literature concerning basis-invariant 
signals of discrete symmetries in the THDM. Thus, in the ERPS = is a necessary and sufficient condition for the 
presence of a U{1) symmetry. 

III. VACUUM STRUCTURE AND RENORMALIZATION 

The presence of a U{1) symmetry in the Higgs potential may (or not) imply the existence of a massless scalar, 
the axion, depending on whether (or not) the U{1) is broken by the vevs. In the previous section we related the 
basis-invariant condition D — in the ERPS with the presence of a U{1) symmetry. In this section we will show that, 
whenever the basis-invariant condition D = is satisfied in the ERPS, there is always a stationary point for which a 
massless scalar, other than the usual Goldstone bosons, exists. 

We start by writing the extremum conditions for the THDM in the ERPS. For simplicity, we will be working in a 
basis where all the parameters are real Q. From Eqs. ([5]) and ([8]), we obtain 

0= ^111^1 + ^ [Xivf + A345W1W2 + XG{3vfv2 - W2)] , 

= Yii W2 + ^ [Ai vl + A345 V2 v\ + Ae (y\ - ivl wi)] , (48) 

where we have defined A345 = A3 -I- A4 -|- A5. We now compute the mass matrices. As we will be considering only 
vacua with real vevs, there will be no mixing between the real and imaginary parts of the doublets. As such, we can 
define the mass matrix of the CP-even scalars as given by 

\Ml\ = - ^— f49) 

L ''J^J 2 9Re($f') aRe($0) ^ ' 

where is the neutral (lower) component of the <I>i doublet. Thus, we obtain, for the entries of this matrix, the 
following expressions: 

[M^]^^ = Yii + i (3Aiw? -f A345 i'2 + 6 Ae 1^1 1^2) 



[M^]22 =yii + 2 (3Aiw2 + A345'yi + 6 Ae 1^2) 

[Af2] = A345 v^V2 Ae {v\ ~vl) . (50) 
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Likewise, the pseudoscalar mass matrix is defined as 



whose entries are given by 



1 d^V 

My = 2 aim($0)aim($0) ^^^^ 

[Ml]^^ =Yn + ^ [Xivl + (A3 + A4 - A5) + 2XeViV2] 
[Ml] 22 ^Yn + ^ [Xivl + (A3 + A4 - As) vf - 2 Ag wi va] 

[Ml] ^2 =>^5ViV2 + ^Xeivl -vl) . (52) 
The expressions (|50p and ([5^ are vahd for all the particular cases we will now consider. 

A. Case = 0, {vi , V2} 7^ 

Let us first study the case Ag = 0, wherein we may solve the extremum conditions in an analytical manner. It is 
trivial to see that Eqs. (|48|) have three types of solutions: both vevs different from zero, one vev equal to zero (say, 
V2) and both vevs zero (trivial non-interesting solution). For a solution with {vi , V2} 7^ 0, a necessary condition 
must be obeyed so that there is a solution to Eqs. 



Ai ~ Xl,, ^ . (53) 
If we use the extremum conditions to evaluate [M|] , we obtain 

[Mi] = f/^"' A345 .l.2\ (54) 

L \X345viv2 Xivi J ^ ^ 

which only has a zero eigenvalue if Eq. (|53p is broken. Thus, there is no axion in this matrix in this case. As for 
[Ml] , we get 

which clearly has a zero eigenvalue corresponding to the Z Goldstone boson. Further, this matrix will have an axion 
if A5 = 0, which is the first condition of Eq. (jiTj) . 



B. Case Ae = 0, {vi / 0, U2 = 0} 
Returning to Eq. ([35]), this case gives us 



Yii = ~\ Ai vf , (56) 



which implies Yn < 0. With this condition, the mass matrices become considerably simpler: 

ivf 

i(A345 -Ai)w2 



[^^]-Ko (Aa+A^As-A.).? ■ ^''^ 



and 



So, we can have an axion in the matrix (j57p if 

A345 - Ai = ^ As = Ai - A3 - A4 (59) 

or an axion in matrix (j58p if 

As = -Ai + A3 + A4 . (60) 
That is, we have an axion if the second or third conditions of Eq. (|Tf|) are satisfied. The other possible case, 



{vi = 0, V2 7^ 0} , produces exactly the same conclusions. 
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C. Case Ae / 

This is the hardest case to treat, since we cannot obtain analytical expressions for the vevs. Nevertheless a full 
analytical treatment is still possible. First, notice that with Ag 7^ Eqs. (HS)) imply that both vevs have to be non- 
zero. At the stationary point of Eqs. the pseudoscalar mass matrix has a Goldstone boson and an eigenvalue 
given by 

-X,[vl+vl)-X,i—^ . (61) 

2. Vi V2 



So, an axion exists if we have 



vf - vl ^ 2 As 
vi V2 Ae 

On the other hand, after some algebraic manipulation, it is simple to obtain from (|48p the following condition: 



(62) 



Ai — A345 — Ag — o (po) 



Substituting Eq. (IS^ into we obtain 



Ai - A345 = ^6 ^ + ^ j ^ 2A2 = A5(Ai - A3 - A4 + A5). (64) 

Thus, we have shown that all of the conditions stemming from the basis-invariant condition D — Q guarantee the 
existence of some stationary point for which the scalar potential yields an axion. Notice that, however, this stationary 
point need not coincide with the global minimum of the potential. 

D. Renormalization group invariance 

We now briefly examine the renormalization group (RG) behavior of our basis-invariant condition D = Q. It would 
be meaningless to say that D — Q implies a J7(l) symmetry if that condition were only valid at a given renormalization 
scale. That is, it could well be that a numerical accident forces D — Q aX only a given scale. To avoid such a conclusion, 
we must verify if D = is a RG-invariant condition (in addition to being basis-invariant). For a given renormalization 
scale /i, the /3- function of a given parameter x is defined as (3^ — ^Ji. dx/d^i. For simplicity, let us rewrite D in Eq. (j45p 
as 

D = -^DlD2 , (65) 



with 



Di — A5(Ai — A3 — A4 + A5) — 2Ag 

D2= (Ai - A3 - A4 - As)" + 16A2 . (66) 



If we apply the operator ^d/d^ to D, we obtain 

Pd = - ^{2D^ D2 Pd^ + Dl Pd,) ■ (67) 

If Z?! = (which corresponds to three of the conditions presented in Eqs. and (HZ])) then we immediately have 
(3d = 0. That is, if D = at a given scale, it is zero at all scales. 

If D2 = and Dl 7^ we will only have /3d — if I3d2 = 0, or equivalently, 

2 (Ai - A3 - A4 - As) (/3ai - /3a3 - Px, - Px,) + 32 Pxe Ag = • (68) 

Given that D2 = implies that Ag = and As = Ai — A3 — A4, we once again obtain Pd = 0. 

Thus, the condition D = is RG-invariant. A direct verification of the RG invariance of Eqs. and (ITT)) , and of 
the conditions that define the ERPS itself, would require the explicit form of the /3 functions of the THDM involving 
the Ag coupling. That verification will be made elsewhere [26j . 
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IV. GENERALIZED CP SYMMETRIES 

It is common to consider the standard CP transformation of the scalar fields as 

^a{t,x)~^^^''{t,x)^<p:{t,~x), (69) 

where the reference to the time (t) and space (x) coordinates will henceforth be suppressed. However, in the presence 
of several scalars with the same quantum numbers, basis transformations can be included in the definition of the CP 
transformation. This yields generalized CP transformations (GCP), 

= x:^'s>i^x:j^ir, (70) 

where X is an arbitrary unitary matrix p?!. [28l|.^ 

Note that the transformation <i>a where is given by Eq. (|70p . leaves the kinetic terms invariant. 

The GCP transformation of a field bilinear yields 

$tGCp$?cp = x:,x,,($„$t)T, (71) 

Under this GCP transformation, the quadratic terms of the potential may be written as 

n^^tGCP^GCP ^ YatX:^X,p<^>l<^>^ 

= x,pY,ix:^^l^^ 

- X^aY:^X;,^i^, = (Xt Y X):,^l^,. (72) 

We have used the Hermiticity condition Yat = i^^* in going to the second line; and changed the dummy indices a ^ (3 
and 6 <-> a in going to the third line. A similar argument can be made for the quartic terms. We conclude that the 
potential is invariant under the GCP transformation of Eq. (|70)) if and only if the coefficients obey 



^ab — Xl^^YapXpb — [X'^ Y X)ab, 
Kb,cd = X^^^X*^Zap^^5XpbXsd- (73) 

Introducing 

AYab = Yab - X^aY^pXlb = [Y ~ {X^ Y Xf] , 
AZab.cd — Zab,cd ^ XaaXjcZ^p^^gXp^Xg^^. (74) 

we may write the conditions for invariance under GCP as 

AYab = 0, (75) 
AZab,cd = 0. (76) 



Given Eqs. (|4|), it is easy to show that 



AYab = An;, 

AZab.cd = AZcd,ab — AZ^a.dc- (^7) 



Thus, we need only consider the real coefficients AFn, AY22, AZu^n, AZ22,22, AZii_22, AZi2,2i, and the complex 
coefficients AY12, AZii^i2, AZ22,i2, and AZi2,i2. 



Equivalently, one can consider a generalized time-reversal transformation proposed in Ref. [29j and considered further in Appendix A 
of Ref. [|. 
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A. GCP and basis transformations 

We now turn to the interplay between GCP transformations and basis transformations. Consider the potential of 
Eq. ^ and call it V{^). Now consider the potential obtained from V{^) by the basis transformation ^'a ~ 

= y:,{^'^^',) + ^z'^,,,A<K)i^':'^',), (78) 

where the coefficients in the new basis are given by Eqs. (1161) and (|17p . We will now prove the following theorem: If 
V{^) is invariant under the GCP transformation of Eq. (|70|l with the matrix X, then V{^') is invariant under a new 
GCP transformation with matrix 

X' = UXU^. (79) 

By hypothesis V{^) is invariant under the GCP transformation of Eq. ((70l) with the matrix X. Eq. ((73)) guarantees 
that Y* = X'^YX. Now, Eq. (fT6|) relates the coefficients in the two basis through Y = U'^Y'U. Substituting gives 

U^Y'*U* = Xt([/ty'J7)X, (80) 

or 

Y'* = {U*X''U'')Y'{UXU^) = X'^Y'X', (81) 

as required. A similar argument holds for the quartic terms and the proof is complete. 

The fact that the transpose C/^ appears in Eq. ([79]) rather than W is crucial. In Eq. l[23|) . applicable to HF 
symmetries, W appears. Consequently, a basis may be chosen where the HF symmetry is represented by a diagonal 
matrix S. The presence of in Eq. (|79p implies that, contrary to popular belief, it is not possible to reduce all GCP 
transformations to the standard CP transformation of Eq. (I69p by a basis transformation. What is possible, as we 
shall see below, is to reduce an invariance of the THDM potential under any GCP transformation, to an invariance 
under the standard CP transformation plus some extra constraints. 

To be more specific, the following result is easily established. If the unitary matrix X is symmetric, then it follows 
that^ a unitary matrix U exists such that X' — UXU^ = 1, in which case Y'* = Y' . In this case, a basis exists 
in which the GCP is a standard CP transformation. In contrast, if the unitary matrix X is not symmetric, then no 
basis exists in which Y and Z arc real for generic values of the scalar potential parameters. Nevertheless, as we shall 
demonstrate below, by imposing the GCP symmetry on the scalar potential, the parameters of the scalar potential 
are constrained in such a way that for an appropriately chosen basis change, Y'* = X''^Y'X' = Y' (with a similar 
result for Z'). 

GCP transformations were studied in Refs. [13, H^. In particular, Ecker, Grimus, and Neufeld [1^ proved that for 
every matrix X there exists a unitary matrix U such that X' can be reduced to the form 

X'^UXU^=( ^^"JV (82) 

y — sm 6^ cos J ' ^ ' 

where < < 7r/2. Notice the restricted range for 0. The value of 9 can be determined in either of two ways: (i) the 
eigenvalues of {X + X"^)'^{X + X^)/2 are cos0, each of which is twice degenerate; or (ii) XX* has the eigenvalues 



B. The three classes of GCP symmetries 

Having reached the special form of X' in Eq. ([5^ . we will now follow the strategy adopted by Ferreira and Silva 
[23} in connection with HF symmetries. We substitute Eq. ([82)) for X in Eq. ([73)) . in order to identify the constraints 
imposed by this reduced form of the GCP transformations on the quadratic and quartic couplings. For each value of 
9, certain constraints will be forced upon the couplings. If two different values of 9 enforce the same constraints, we 
will say that they are in the same class (since no experimental distinction between the two will then be possible) . We 
will start by considering the special cases of 6* = and 9 = 7r/2, and then turn our attention to < < 7r/2. 



Here, we make use of a theorem in linear algebra that states that for any unitary symmetric matrix X, a unitary matrix V exists such 
that X = VV^ . A proof of this result can be found, e.g., in Appendix B of Ref. Q. 
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1. CPl: 6^0 



When — 0, X' is the unit matrix, and we obtain the standard CP transformation, 



$1 ^ $1, 

$2 ^ (83) 



under which Eqs. (|73p take the very simple form 



V* — V , 

^ ab ~ ^abi 

Kb.cd = Zab^cd- (84) 

We denote this CP transformation by CPl. It forces aU couphngs to be real. Since most couplings are real by the 
Hermiticity of the Higgs potential, the only relevant constraints are Im 7/1^2 = ImAs = ImAg = ImAr = 0. 



2. CP2: e = Tv/2 



When = tt/2, 



and we obtain the CP transformation, 



X' 



(85) 



$1 

$2 



^2, 



(86) 



which we denote by CP2. This was considered by Davidson and Haber [s*] in their Eq. (37), who noted that if this 
symmetry holds in one basis, it holds in all basis choices. Under this transformation, Eq. (|75p forces the matrix of 
quadratic couplings to obey 



= Ar 



- m2 
2* 
2 



-2vn\^ 

9 2 



(87) 



leading to m22 — and rn\^ = 0. Similarly, we may construct a matrix of matrices containing all coefficients 
AZab^cd- The uppermost-leftmost matrix corresponds to AZu^cd- The next matrix along the same line corresponds 
to AZi2,c(i) and so on. To enforce invariance under CP2, we equate it to zero. 



= 



/ / Ai — A2 Ag + A7 
A^ + A? 



A^ + A^ 
A^ + A? 



Ae + Ay 






Ae + At 



Ae + A7 

Ag + Ay A2 — Ai y y 



We learn that invariance under CP2 forces rtv^^ — "Jii and rn\^ = 0, A2 = Ai, and Ay = — Ag, leading precisely to the 
ERPS of Eq. ([55)1 . Recall that Gunion and Haber Q found that, under these conditions we can always find a basis 
where all parameters are real. As a result, if the potential is invariant under CP2, there is a basis where CP2 still 
holds and in which the potential is also invariant under CPl. 



3. CPS: < 6* < 7r/2 

Finally we turn to the cases where < 6* < 7r/2. Imposing Eq. ([75]) yields 

= AYii = [(m^j - 77122) s - 2 Rem^2 c] s, 
= Ay22 - -AYii, 

= AY12 = Re 777 i2 (c2 - 1) - 27 Im 777 i2 + 5(^22 " "^\\) ^s, (89) 
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where we have used c = cos 6*, s = sin 61, C2 = cos 20, and S2 = sin 20. Since 6 ^ 0, 7r/2, the conditions TO22 
= are imposed, as in CP2. Similarly, Eq. (|76p yields 



and 



= AZnai 

= AZ22,22 

- AZn,22 

= AZi2,21 

= RcAZn.i2 

= ReAZ22,i2 

= RcAZi242 
= Im AZii,i2 

~ Im AZ22,i2 
= Im AZ12.12 



= Ai(l - c") - A2s'^ - iA345S2 +4 ReA6C^s + 4 ReAycs^, 
^ A2(l - c^) - Xis^ - 5A345S2 - 4 Re Ayc^s - 4 ReA6Cs^ 



-is2 [4Re (As - A7)c2 + (Ai + A2 
AZii 22 

is [(-3Ai + A2 + 2A345)c - (Ai + A2 



+4Re A6(2s 
is[(-Ai+3A2- 

-4ReA6S3 4 
AZii 22 

i[ImA6(3 + C2)- 
i[ImA6(l-C2)- 



f S3) -4ReA7S3] , 

2A345)C+ (Ai + A2 

4ReA7(2s + S3)], 



2A345)S2] , 
2A345)C3 
2A345)C3 



= 2c[ImAr, 



ImA7(l - C2) - ImA5S2] , 
Im A7(3 + C2) + Im A5S2] , 
Im(A6 - A7)s] , 



where A345 = A3 + A4 + Re A5, C3 = cos 30, and S3 = sin 30. 
The last three equations may be written as 



= 



(1 



'S2 
S2 
f C2) 



(3 + C2) 
(I-C2) 

S2 



(1 

(3- 



~C2) 

I-C2) 

-32 



Im A5 
Im Xq 
Im A7 



(90) 



(91) 



(92) 



The determinant of this homogeneous system of three equations in three unknowns is 32c^, which can never be zero 
since we are assuming that ^ 7r/2. As a result, A5, Ag, and A7 are real, whatever the value ofO<0<7r/2 chosen 
for the GCP transformation. Since to^2 — 0? potentially complex parameters must be real. We conclude that a 
potential invariant under any GCP with O<0<7r/2is automatically invariant under CPl. Combining this with what 
we learned from CP2, we conclude the following: if a potential is invariant under some GCP transformation, then a 
basis may be found in which it is also invariant under the standard CP transformation, with some added constraints 
on the parameters. 

The other set of five independent homogeneous equations in five unknowns has a determinant equal to zero, meaning 
that not all parameters must vanish. We find that 



— AZii 11 — AZ22,22 

= ReAZii42-ReAZ22,i2 
Since s 7^ 0, we obtain the homogeneous system 



2s [s (Ai - 
s [-C (Ai 



A2) 
- A2) 



c 2Re (Ag - 
- s 2Re (Ae 



Ar)], 
-At)]. 







s 

— c 



Ai — A2 
2Re (Ae + A7) 



(93) 



(94) 



Thus, GCP invariance with any value of 
AZ22.22 = -AZii 92 and 



whose determinant is unity. We conclude that A2 = Ai and A7 = — Ag 
< < tt/2 leads to the ERPS of Eq. Substituting back we obtain AZn^n 

Re AZ11.12 = — Re AZ22,i2, leaving only two independent equations: 

= AZii^ii = is2 [(Ai - A345)S2 + 4A6C2] , 
= ReAZ22,12 = ^S2 [(Ai - A345)C2 - 4A6S2] , (95) 

where we have used c + C3 = 2cc2 and s + S3 = 2cs2. Since S2 7^ 0, the determinant of the system does not vanish, 
forcing Ai = A345 and Ag = 0. 

Notice that our results do not depend on which exact value ofO<0<Tr/2in Eq. ([82| we have chosen. If we require 
invariance of the potential under GCP with some particular value of < < 7r/2, then the potential is immediately 
invariant under GCP with any other value of < < tt/2. We name this class of CP invariances, CP3. Combining 
everything, we conclude that invariance under CP3 implies 



mil = "^2 

A2 = Ai, A7 = Ag = 0, 

Im A5 = 0, Re A5 = Ai — A3 

The results of this section are all summarized in Table [T] of section |Vl 



^22 1 



mi2 = 0, 



A4. 



(96) 
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C. The square of the GCP transformation 

If we apply a GCP transformation twice to the scalar fields, we will have, from Eq. ([70]). that 

(^GCP)GCP ^ x^^^^Gcy ^ x^^x:,t>, , (97) 
so that the square of a GCP transformation is given by 

{GCPf ^ XX* . (98) 

In particular, for a generic unitary matrix X , (GCP)^ is a Higgs Family symmetry transformation. 

Usually, only GCP transformations with (GCP)^ — 1 (where 1 is the unit matrix) are considered in the literature. 
For such a situation, X = = X*, and one can always find a basis in which X = 1. In this case, a GCP 
transformation is equivalent to a standard CP transformation in the latter basis choice. For example, the restriction 
that [GCPY = 1 (or equivalently, requiring the squared of the corresponding generalized time-reversal transformation 
to equal the unit matrix) was imposed in Ref. 9] and more recently in Ref. 15]. However, as we have illustrated in 
this section, the invariance under a GCP transformation, in which (GCP)^ ^ 1 (corresponding to a unitary matrix 
X that is not symmetric) is a stronger restriction on the parameters of the scalar potential than the invariance under 
a standard CP transformation. 

As we see from the results in the previous sections, X is not symmetric for the symmetries CP2 and CPS. In fact, 
this feature provides a strong distinction among the three GCP symmetries previously introduced. Let us briefly 
examine (GCP)^ for the three possible cases CPl, CP2 and CPS. 

1. {CPlf 

Comparing Eqs. ([70)1 and ([55]) . we come to the immediate conclusion that Xcpi = 1, so that Eq. ([M)) yields 

{CPlf = 1. (99) 
This implies that a CPl-invariant scalar potential is invariant under the symmetry group Z2 = {1 , CPl}. 



2. {CP2f 

The matrix Xcp2 is shown in Eq. ((85|) so that, by Eq. (f98|) . we obtain 

(CP2)2 = -1. (100) 

Although this result significantly distinguished CP2 from CPl, the authors of Ref. J^] noted (in considering their 

CPg^^ symmetries) that the transformation law for under (CP2)^ can be reduced to the identity by a global 
hypercharge transformation. That is, if we start with the symmetry group Z4 — {1 , CP2 , —1, — CP2}, we can 
impose an equivalence relation by identifying two elements of related by multiplication by —1. If we denote 
{Z2)y = {1 , —1} as the two-element discrete subgroup of the global hypercharge U{1)y, then the discrete symmetry 
group that is orthogonal to U{1)y is given by ^4/(^2)1- = ^2- Hence, the CP2-invariant scalar potential exhibits a 
Z2 symmetry orthogonal to the Higgs flavor symmetries of the potential. 



3. [CPif 

The matrix ATcps is given in Eq. ([5^ . with < < 7r/2, so that, by Eq. we obtain 

which once again is not the unit matrix. However, the transformation law for under (CPS)^ cannot be reduced to 
the identity by a global hypercharge transformation. This is the reason why Ref. [15] did not consider CPS. However, 
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(CP3)^ is a non-trivial HF symmetry of the CP3-invariant scalar potential.^ Thus, one can always reduce the square 
of CPS to the identity by applying a suitable HF symmetry transformation. In particular, a CP3-invariant scalar 
potential also exhibits a Zi symmetry that is orthogonal to the Higgs flavor symmetries of the potential. 

In this paper, we prove that there are three and only three classes of GCP transformations. Of course, within each 
class, one may change the explicit form of the scalar potential by a suitable basis transformation; but that will not 
alter its physical consequences. Similarly, one can set some parameters to zero in some ad-hoc fashion, not rooted in 
a symmetry requirement. But, as we have shown, the constraints imposed on the scalar potential by a single GCP 
symmetry can be grouped into three classes: CPl, CP2, and CPS. 

V. CLASSIFICATION OF THE HF AND GCP TRANSFORMATION CLASSES IN THE THDM 

A. Constraints on scalar potential parameters 

Suppose that one is allowed one single symmetry requirement for the potential in the THDM. One can choose an 
invariance under one particular Higgs Family symmetry. We know that there are only two independent classes of such 
simple symmetries: Zi and Peccei-Quinn J7(l). One can also choose an invariance under a particular GCP symmetry. 
We have proved that there are three classes of GCP symmetries, named CPl, CP2, and CPS. If any of the above 
symmetries is imposed on the THDM scalar potential (in a specified basis), then the coefficients of the scalar potential 
are constrained, as summarized in Table U For completeness, we also exhibit the constraints imposed by 50(3), the 
largest possible continuous HF symmetry that is orthogonal to the global hypercharge [/(l)y transformation. 





TABLE I: Impact of the symmetries 


on the coefficients of the Higf 


;s potential in a specified basis. 






symmetry 




mi2 Ai 


A2 A3 A4 


As 


As 


At 


Z2 

















17(1) 


















50(3) 







Ai Ai — A3 











Ha 


mfi 


real 


Ai 


real 




Ae 


CPl 




real 




real real 


real 


CP2 


mfi 





Ai 






-As 


CPS 







Ai 


Ai — A3 — A4 (real) 









Empty entries in Table U correspond to a lack of constraints on the corresponding parameters. Table [TJias been 
constructed for those basis choices in which Z2 and C/(l) have the specific forms in Eqs. ([24]) and (|27|) . respectively. 
If, for example, the basis is changed and Z2 acquires the form H2 in Eqs. (P5|) . then the constraints on the coefficients 
are altered, as shown explicitly on the fourth line of Table HI However, this does not correspond to a new model. All 
physical predictions are the same since the specific forms of Zi and H2 differ only by the basis change in Eq. (|26p . 
The constraints for CPl, CP2, and CPS shown in Table I apply to the basis in which the GCP transformation of 
Eq. ([70]) is used where X has been transformed into X' given by Eq. ([82]) . with = 0, = 7r/2, and < < 7r/2, 
respectively. 

B. Multiple symmetries and GCP 

We now wish to consider the possibility of simultaneously imposing more than one symmetry requirement on the 
Higgs potential. For example, one can require that Zi and H2 be enforced within the same basis. In what follows, we 
shall indicate that the two symmetries are enforced simultaneously by writing © H2. Combining the constraints 
from the appropriate rows of Table IH we conclude that, under these two simultaneous requirements 

22 2 n 

A2 = Ai, A7 = A6=0, ImA5==0. (102) 



In Section Ivb. we shall identify (CP3)^ with the Peccei Quinn U(l) symmetry defined as in Eq. I|27^ and then transformed to a new 
basis according to the unitary matrix defined in Eq. 111051 . 
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This coincides exactly with the conditions of the ERPS in a very special basis, as shown in Eq. p4|) . Since CP2 leads 
to the ERPS of Eq. we conclude that 

^2 © 112 = CP2 in some specific basis. (103) 

This was noted previously by Davidson and Haber ^8(| . Now that we know what all classes of HE and CP symmetries 
can look like, we can ask whether all GCP symmetries can be written as the result of some multiple HF symmetry. 

This is clearly not possible for CPl because of parameter counting. Table |T] shows that CPl reduces the scalar 
potential to ten real parameters. We can still perform an orthogonal basis change while keeping all parameters real. 
This freedom can be used to remove one further parameter; for example, setting 771^2 = by diagonalizing the Y 
matrix. No further simplification is allowed. As a result, CPl leaves nine independent parameters. The smallest HE 
symmetry is Z2. Table H] shows that Z2 reduces the potential to six real and one complex parameter. The resulting 
eight parameters could never account for the nine needed to fully describe the most general model with the standard 
CP invariance CPl.^ 

But one can utilize two HE symmetries in order to obtain the same constraints obtained by invariance under CPS. 
Let us impose both U{1) and n2 in the same basis. From Table [H we conclude that, under these two simultaneous 
requirements 

22 2 n 

™22=™ll7 'Tlia = 0, 

A2 = Ai, A7 = A6 = 0, A5=0. (104) 

This does not coincide with the conditions for invariance under CPS shown in Eq. (|96p . However, one can use the 
transformation rules in Eqs. (A1S)-(A2S) of Davidson and Haber in order to show that a basis transformation, 

may be chosen which takes us from Eqs. ([96|) . where ReAs ~ Ai — A3 — A4, to Eqs. (|104p . where A5 = (while 
maintaining the other relations among the scalar potential parameters). We conclude that 

U{1) H2 = CPS in some specific basis. (106) 

Note that in the basis in which the CPS relations of Eq. (|96p are satisfied with A5 ^ 0, the discrete HF symmetry 
H2 is still respected. However, using Eq. ()105|) . it follows that the U(l)-Peccei Quinn symmetry corresponds to the 
invariance of the scalar potential under <1>q Oab^b, where O is an arbitrary 5*0(2) matrix. 

The above results suggest that it should be possible to distinguish CPl, CP2, and CPS in a basis invariant fashion. 
Botella and Silva Q have built three so-called J-invariants that detect any signal of CP violation (either explicit 
or spontaneous) after the minimization of the scalar potential. However, in this paper we are concerned about the 
symmetries of the scalar potential independently of the choice of vacuum. Thus, we shall consider the four so-called 
J-invariants built by Gunion and Haber [§] in order to detect any signal of explicit CP violation present (before the 
vacuum state is determined). If any of these invariants is nonzero, then CP is explicitly violated, and neither CPl, 
nor CP2, nor CPS hold. Conversely, if all /-invariants are zero, then CP is explicitly conserved, but we cannot tell a 
priori which GCP applies. Eqs. (|10Sp and p06p provide the crucial hint. If we have CP conservation, Z2 © H2 holds, 
and U{1) does not, then we have CP2. Alternatively, if we have CP conservation, and U{1) © H2 also holds, then we 
have CPS. We recall that both CP2 and CPS lead to the ERPS, and that the general conditions for the ERPS in 
Eq. are basis independent. This allows us to distinguish CP2 and CPS from CPl. But, prior to the present work, 
no basis-independent quantity had been identified in the literature that could distinguish Z2 and C/(l) in the ERPS. 
The basis-independent quantity D introduced in subsection III El is precisely the invariant required for this task. That 
is, in the ERPS D^O imphes CP2, whereas D = imphes CPS. 

One further consequence of the results of Table U can be seen by simultaneously imposing the U(l) Peccei-Quinn 
symmetry and the CPS symmetry in the same basis. The resulting constraints on the scalar potential parameters are 
precisely those of the SO(S) HF symmetry. Thus, we conclude that 

[/(I) © CPS = 50(3). (107) 

In particular, SO{S) is not a simple HF symmetry, as the invariance of the scalar potential under a single element of 
SO(S) is not sufficient to guarantee invariance under the full SO(S) group of transformations. 



In Ivanov's language, this is clear since CPl corresponds to a Z2 transformation of the vector r, which is the simplest transformation 
on rone could possibly make. See section [VDI 
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C. Maximal symmetry group of the scalar potential orthogonal to U{1)y 

The standard CP symmetry, CPl, is a discrete Z2 symmetry that transforms the scalar fields into their complex 
conjugates, and hence is not a subgroup of the U{2) transformation group of Eq. [151 We have previously noted 
that THDM scalar potentials that exhibit any non-trivial HF symmetry G is automatically CP-conserving. Thus, 
the actual symmetry group of the scalar potential is in fact the semidirect product* of G and Z2, which we write 
as G X Z2. Noting that U{1) x Z2 = 50(2) x = 0(2), and 5*0(3) Z2 = 0(3), we conclude that the maximal 
symmetry groups of the scalar potential orthogonal to U{1)y for the possible choices of HF symmetries are given in 
Table 



TABLE 11: Maximal symmetry groups [orthogonal to global U{1)y hypercharge] of the scalar sector of the THDM. 
designation HF symmetry group maximal symmetry group 



Z2 Z2 iJZ; 



^2 



2J 

Peccei-Quinn (7(1) 0(2) 

5*0(3) SO(?,) 0(3) 

CPl — Z2 

CP2 {Z^f {Z-if 

CPS 0(2) 0(2) (g) Z-2 



Finally, we reconsider CP2 and CPS. Eq. (|103|) implies that the CP2 symmetry is equivalent to a {Z^f' HF 
symmetry. To prove this statement, we note that in the two-dimensional flavor space of Higgs fields, the Z2 and n2 
discrete symmetries defined by Eqs. (|24p and (PS)) are given by: 

Z2 = {So , 5i} , n2 = {So , 52} , (108) 

where 5o = 1 is the 2x2 identity matrix and 




(109) 



If we impose the Z2 and n2 symmetry in the same basis, then the scalar potential is invariant under the dihedral 
group of eight elements, 

D4 = {So , Si , S2 , S3 , —So , —Si , —52 , —53} , (110) 

where 53 = 5i52 = — 525i. As before, we identify {Z2)y = {So, —So} as the two-element discrete subgroup of 
the global hypercharge J7(l)y. However, we have defined the HF symmetries to be orthogonal to U{1)y- Thus, 
to determine the HF symmetry group of CP2, we identify as equivalent those elements of that are related by 
multiplication by —So- Group theoretically, we identify the HF symmetry group of CP2 as 

Di/{Z2)Y ^Z2®Z2. (Ill) 

The HF symmetry group of CP2 is not the maximally allowed symmetry group. In particular, the constraints of 
CP2 on the scalar potential imply the existence of a basis in which all scalar potential parameters are real. Thus, the 
scalar potential is explicitly CP-conserving. The Z2 symmetry associated with this CP transformation is orthogonal 
to the HF symmetry as previously noted. (This is easily checked explicitly by employing a four-dimensional real 
representation of the two complex scalar fields.) Thus, the maximal symmetry group of the CP2-symmetric scalar 
potential is {Z2)^ . Similarly, Eq. (I106|) implies that the CP3 symmetry is equivalent to a IJ(\) xi Z2 HF symmetry. 
This is isomorphic to an 0(2) HF symmetry, which is a subgroup of the maximally allowed 50(3) HF symmetry 
group. However, the constraints of CP3 on the scalar potential imply the existence of a basis in which all scalar 



In general, the non-trivial element of Z2 will not commute with all elements of G, in which case the relevant mathematical structure is 
that of a semidirect product. In cases where the non-trivial element of Z2 commutes with all elements of G, we denote the corresponding 
direct product as G ^ ^2 • 

For ease of notation, we denote Z2 ® Z2 by (.^2)^ and Z2 ® Z2 ® Z2 by {Z2)^- 
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potential parameters are real. Thus, the scalar potential is explicitly CP-conserving. Once again, the Zi symmetry 
associated with this CP transformation is orthogonal to the HP symmetry noted above. Thus, the maximal symmetry 
group of the CP3-symmetric scalar potential is 0(2) ® Z2. 

The above results are also summarized in Table HIl In all cases, the maximal symmetry group is a direct product of 
the HP symmetry group and the Zi corresponding to the standard CP-transformation, whose square is the identity 
operator. 

One may now ask whether Table [IT] exhausts all possible independent symmetry constraints that one may place 
on the Higgs potential. Perhaps one can choose other combinations, or maybe one can combine three, four, or more 
symmetries. We know of no way to answer this problem based only on the transformations of the scalar fields $a- 
Fortunately, Ivanov has solved this problem Q by looking at the transformation properties of field bilinears, thus 
obtaining for the first time the list of symmetries given in the last column of Table |TT1 



D. More on multiple symmetries 



We start by looking at the implications of the symmetries we have studied so far on the vector f = {t"!, ^2, rs}, 
whose components were introduced in Eq. (fTO|) . Notice that a unitary transformation U on the fields $a induces an 
orthogonal transformation O on the vector of bilinears r, given by Eq. (|37p . For every pair of unitary transformations 
±J7 of SU(2), one can find some corresponding transformation O of S0{3), in a two-to-one correspondence. We then 
see what these symmetries imply for the coefficients of Eq. ([SJ (recall the A^j, is a symmetric matrix). Below, we list 
the transformation of r under which the scalar potential is invariant, followed by the corresponding constraints on 
the quadratic and quartic scalar potential parameters, and A^^. 

Using the results of Table I, we find that Z2 implies 



-n 



U{1) implies 
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(112) 



(113) 



(114) 



where O is an arbitrary 3x3 orthogonal matrix of unit determinant. In the language of bilinears, a basis invariant 
condition for the presence of iS'0(3) is that the three eigenvalues of A are equal. (Recall that A = {A^j}; i,j = 1, 2, 3). 
As for the GCP symmetries, CPl implies 



CP2 implies 
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(115) 



(116) 
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and CPS implies 
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r — > 
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(117) 



Notice that in CP3 two of the eigenvalues of A are equal, in accordance with our observation that D can be used to 
distinguish between CP2 and CP3. 

Because each unitary transformation on the fields $a induces an S0{3) transformation on the vector of bilinears r, 
and because the standard CP transformation corresponds to an inversion of r2 (a Z2 transformation on the vector r), 
Ivanov considers all possible proper and improper transformations of 0(3) acting on r. He identifies the following 
six classes of transformations: (i) Z2; (ii) (^2)^; (iu) (-^2)^; (iv) 0(2); (v) 0(2) (g) Z2; and (vi) 0(3). Note that these 
symmetries are all orthogonal to the global U{1)y hypercharge symmetry, as the bilinears and r are all singlets 
under a U{1)y transformation. The six classes above identified by Ivanov correspond precisely to the six possible 
maximal symmetry groups identified in Table |TT1 No other independent symmetry transformations are possible. 

Our work permits one to identify the abstract transformation of field bilinears utilized by Ivanov in terms of 
transformations on the scalar fields themselves, as needed for model building. Combining our work with Ivanov's, we 
conclude that there is only one new type of symmetry requirement which one can place on the Higgs potential via 
multiple symmetries. Combining this with our earlier results, we conclude that all possible symmetries on the scalar 
sector of the THDM can be reduced to multiple HF symmetries, with the exception of the standard CP transformation 
(CPl). 



VI. BUILDING ALL SYMMETRIES WITH THE STANDARD CP 

We have seen that there are only six independent symmetry requirements, listed in Table |TT1 that one can impose 
on the Higgs potential. We have shown that all possible symmetries of the scalar sector of the THDM can be reduced 
to multiple HF symmetries, with the exception of the standard CP transformation (CPl). Now we wish to show a 
dramatic result: all possible symmetries on the scalar sector of the THDM can be reduced to multiple applications of 
the standard CP symmetry. 

Using Eq. ((79|) , we see that the basis transformation of Eq. (fT5|) , changes the standard CP symmetry of Eq. ((69|) 
into the GCP symmetry of Eq. ^0^, with 

X ^ UU^. (118) 

In particular, an orthogonal basis transformation does not affect the form of the standard CP transformation. Since 
we wish to generate X ^ 1, we will need complex matrices U . 

Now we wish to consider the following situation. We have a basis (call it the original basis) and impose the standard 
CP symmetry CPl on that original basis. Next we consider the same model in a different basis (call it M) and impose 
the standard CP symmetry on that basis M . In general, this procedure of imposing the standard CP symmetry in 
the original basis and also in the rotated basis M leads to two independent impositions. The first imposition makes 
all parameters real in the original basis. One way to combine the second imposition with the first is to consider the 
basis transformation Um taking us from basis M into the original basis. As we have seen, the standard CP symmetry 
in basis M turns, when written in the original basis, into a symmetry under 

<l>tcP = (XM)L($i)*, (119) 

with Xm = UmUJi- Next we consider several such possibilities. 
We start with 

Ua={ '^/^ -*^V4^ x^=f (120) 

Here and henceforth c (s) with a subindex indicates the cosine (sine) of the angle given in the subindex. We denote 
by CPlyi the imposition of the CP symmetry in Eq. (|119p with Xm = Xa (which coincides with the imposition of 
the standard CP symmetry in the basis M = A). 
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Next we consider 

""-[^ /"')■ Ij^ <'^" 

We denote by CPIb the imposition of the CP symmetry in Eq. (|119p with Xm = Xb (which coincides with the 
imposition of the standard CP symmetry in the basis M ~ B) . 
A third possible choice is 

where 5 ^ n7r/2 with n integer. We denote by CPlc the imposition of the CP symmetry in Eq. pi9p with Xm = Xc 
(which coincides with the imposition of the standard CP symmetry in the basis M = C). 
FinaUy, we consider 

Ud = { V Xu={ (123) 

\ Cs/2 J \ ISS CS J 

where 8 ^ ?i7r/2 with n integer. We denote by GPlo the imposition of the CP symmetry in Eq. (|119p with Xm = Xjj 
(which coincides with the imposition of the standard CP symmetry in the basis M — D). 

The impact of the first three symmetries on the coefficients of the Higgs potential are summarized in Table IIIII 




TABLE III: Impact of the CPIa/ 


symmetries 


on the coefficients of the Hig^ 


;s potential. The notation "imag" means 


that the 


corresponding entry is purely imaginary. CPl 


in the original basis has been 


included for reference. 
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imag 
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imag 


imag 
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lAsle^'* 
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Imposing CP Id on the Higgs potential leads to the more complicated set of equations: 

2Im (m^2) c<5 + (m22 — "^n) S5 = 0, 

2Im (Ae - A7) C25 + A12345 S2S = 0, 
2Im (Ae + Ay) cs + (Ai - A2) S5 0, 

ImAs cs + Re (Ag - A7) sg - 0, (124) 

where 

A12345 = 1 (^1 + ^2) - A3 - A4 + ReAs. (125) 
Combining these results with those in Table [H we have shown that 

CPl ® CPIb = Z2 in some specific basis, 

CPl ©CPlc - C/(l), 

CPl ® CPIa © CPls = CP2 in some specific basis, 

CPl ® CPIa © CPlc = CP3 in some specific basis, 

CPl© CPlc ©CPId = 50(3). (126) 

Let us comment on the "specific basis choices" needed. Imposing CPl © CPIb leads to mfj = Ag = Ay = and 
ImAs — 0, while imposing Z2 leads to m^2 = Ag = Ay = with no restriction on A5. However, when Z2 holds one 
may rephase $2 by the exponential of — i arg(A5)/2, thus making A5 real. In this basis, the restrictions of Z2 coincide 
with the restrictions of CPl © CP Is- Similarly, imposing CPl © CPIa © CPlc leads to 

1X1^2 = and A2 = Ai. We see from Table [I] that CP3 has these features, except that A5 need not vanish; it is real 
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and ReAs = Ai — A3 — A4. Starting from the CPS conditions and using the transformation rules in Eqs. (A13)-(A23) 
of Davidson and Haber ,§3, we find that a basis choice is possible such that ReAs = 0-^° Perhaps it is easier to prove 
the equality 

CPl ® CPIb ® CPlz3 = CP3 in some specific basis. (127) 

In this case, the only difference between the impositions from the two sides of the equality come from the sign of 
ReAs, which is trivial to flip through the basis change $2 —^2- Finally, imposing CPl ® CPIa © CPls we obtain 
= ImAs = Ag = A7 = 0, 777.22 = "^11 ^^'^ A2 = Ai. This does not coincide with the conditions of CP2 which 
lead to the ERPS of Eq. ([55]) . Fortunately, and as we mentioned before, Davidson and Haber 8] proved that one 
may make a further basis transformation such that Eq. (j34p holds, thus coinciding with the conditions imposed by 
CP1©CP1a®CP1b. 

Notice that our description of CP2 in terms of several CPl symmetries is in agreement with the results found by the 
authors of Ref. [lH|. These authors also showed a very interesting result, concerning spontaneous symmetry breaking 
in 2HDM models possessing a CP2 symmetry. Namely, they prove (their Theorem 4) that electroweak symmetry 
breaking will necessarily spontaneously break CP2. However, they also show that the vacuum will respect at least one 
of the CPl symmetries which compose CP2. Which is to say, in a model which has a CP2 symmetry, spontaneous 
symmetry breaking necessarily respect the CPl symmetry. 

In summary, we have proved that all possible symmetries on the scalar sector of the THDM, including Higgs Family 
symmetries, can be reduced to multiple applications of the standard CP symmetry. 

VII. CONCLUSIONS 

We have studied the application of generalized CP symmetries to the THDM, and found that there are only two 
independent classes (CP2 and CPS), in addition to the standard CP symmetry (CPl). These two classes lead to 
an exceptional region of parameter, which exhibits either a Z2 discrete symmetry or a larger C/(l) Peccei-Quinn 
symmetry. We have succeeded in identifying a basis-independent invariant quantity that can distinguish between the 
Z2 and U(l) symmetries. In particular, such an invariant is required in order to distinguish between CP2 and CPS, 
and completes the description of all symmetries in the THDM in terms of basis-invariant quantities. Moreover, CP2 
and CPS can be obtained by combining two Higgs Family symmetries and that this is not possible for CPl. 

We have shown that all symmetries of the THDM previously identified by Ivanov can be achieved through simple 
symmetries, with the exception of SO{3). However, the 5*0(3) Higgs Family symmetry can be achieved by imposing 
a U{1) Peccei-Quinn symmetry and the CPS-symmetry in the same basis. Finally, we have demonstrated that all 
possible symmetries of the scalar sector of the THDM can be reduced to multiple applications of the standard CP 
symmetry. Our complete description of the symmetries on the scalar fields can be combined with symmetries in the 
quark and lepton sectors, to aid in model building. 
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